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ABSTRACT 



We investigated a sample of 15 luminous high-redshift quasars (3.3 < z < 5.1) 
to measure the mass of their super-massive black holes (SMBH) and compare, 
for the first time, results based on Civ, Mgll, and H/? emission lines at high- 
redshifts. Assuming gravitationally bound orbits as dominant broad-line region 
gas motion, we determine black hole masses in the range of M b h ~ 2 x 1O 8 M up 
to M bh ~ 4 x 1O 1O M . While the black hole mass estimates based on Civ and 
H/9 agree well, Mgll typically indicates a factor of ~ 5x lower SMBH masses. A 
flatter slope of the B./3 radius - luminosity relation, a possibly steeper slope of 
the Mg II radius - luminosity relation, and a slightly larger radius of the Mg II 
BLR than for H/5 could relax the discrepancy. In spite of these uncertainties, the 
Civ, Mgll, and H/5 emission lines consistently indicate super-massive black hole 
masses of several times 1O 9 M at redshifts up to z=5.1. Assuming logarithmic 
growth by spherical accretion with a mass to energy conversion efficiency of 
e = 0.1 and an Eddington ratio L bo i/L edd calculated for each quasar individually, 
we estimate black hole growth-times of the order of several ~ 100 Myr which 
are smaller than the age of the universe at the corresponding redshift. Assuming 
high-mass seed black holes (M^ ed = 10 3 tol0 5 M Q ) the SMBHs in the z ~ 3.5 
quasars began to grow at redshifts z > 4, while for the quasars with z > 4.5 
they started at z ~ 6 to 10. These estimated time scales for forming SMBHs at 
high redshifts, together with previous studies indicating high quasar metallicities, 
suggest that the main SMBH growth phase occurs roughly contemporaneously 
with a period of violent and extensive star formation in proto-galactic nuclei. 
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1. Introduction 

During the last few years it has been established that almost all massive galaxies are 
hosting a super-massive black hole (SMBH) in their centers. It has been found that the 
SMHB mass is well correlated with the host galaxy's spheroidal component mass and with 
the bulge luminosity (e.g., Kormendy & Richstone 1995; Magorrian et al. 1998; McLure & 
Dunlop 2001). Furthermore, a remarkably close correlation of the SMBH mass and the stellar 
velocity dispersion, cr*, is detected for active, as well as for non- active galaxies (Gebhardt et 
al. 2000; Graham et al. 2001; Merritt & Ferrarese 2001; McLure k Dunlop 2002; Tremaine et 
al.2002). These results, together with indications that luminous quasars reside in massive 
early type galaxies (e.g., McLure et al. 1999; Kukula et al.2001; Nolan et al. 2001; Percival 
et al.2001; Dunlop et al. 2003), point to a closely related formation of massive spheroidal 
galaxies and the formation and growth of central SMBH. However, it is still under debate 
whether the black hole and the host galaxy formed simultaneously or whether one component 
formed before the other (e.g., Silk & Rees 1998; Cattaneo et al. 1999; Salucci et al. 1999; 
Kauffmann & Haehnelt 2000; Haiman & Loeb 2001; Menou et al. 2001; Archibald et al. 2002; 
DiMatteo et al.2003). 

It is of great importance to estimate SMBH masses at higher redshifts to investigate 
whether the relations of SMBH mass with fundamental parameters of their harboring galaxies 
that are found for active and non-active galaxies in the local universe are already established 
at early cosmic epochs. Recent studies of quasars at moderate (McLure & Jarvis 2002, 
z < 1) and high-redshifts (Shields et al. 2003, z < 3; Vestergaard 2004, z < 6) assume 
that the locally determined radius - luminosity relation can be applied up to redshifts of 
z ~ 6. Particularly, black hole masses have been estimated for several high-redshift quasars, 
providing evidence that even at early cosmic times black holes with several times 10 9 M Q 
are already in place (Barth et al. 2003; Willott et al. 2003; Vestergaard 2004; Warner et 
al.2004). The detection of quasars at redshifts z > 6 (Fan et al.2001) and black hole mass 
estimates well beyond M bh ~ 1O 1O M (Bechthold et al. 2003; Netzer 2003) challenges black 
hole growth models to buildup SMBHs on time scales of several times ~ 100 Myr, because at 
z > 6 the age of the universe is less than 0.9 Gyr. However, the high-redshift Mbh estimates 
are all based on the CivA1549 emission line and the assumption of virialized gas motion. 
It is obviously important to use more than a single emission line to determine Mbh as a 
fundamental AGN parameter. 

Much effort has been spent to derive methods to estimate the black hole mass, Mbh, 
based on reverberation experiments and assuming virial dominated motion of the broad line 
region (BLR) gas (e.g., Kaspi et al. 2000; Peterson & Wandel 1999,2000; Wandel, Peterson, & 
Malkan 1999). This results in an empirically established radius-luminosity relation, R BLR oc 



- 3- 



L@ with e.g. (3 = 0.7 (Kaspi et al. 2000). Since long-term variability studies are currently not 
feasible for quasars, alternative methods are suggested employing single epoch emission line 
profile width measurements to estimate (Vestergaard 2002). There is still some concern 
regarding this approach (e.g., Krolik 2001) and that the R-L relation might be flatter than 
(3 = 0.7 (e.g., Peterson et al. 2000; McLure & Jarvis 2002; Corbett et al. 2003; Netzer 2003; 
Vestergaard 2002, 2004 for a more detailed discussion). However, using emission line profiles 
of H/9 and CivA1549 to estimate the SMBH mass of quasars yields consistent results that 
are on average within ~ 20 % (e.g., Vestergaard 2002; Warner et al. 2003). Recently, McLure 
& Jarvis (2002) suggested MgllA2798 to measure black hole masses. Particularly for quasars 
at redshifts z > 4, Mgll is of great interest because H/3 is shifted out of the K-band. As a low 
ionization line Mg II can be used as a well suited replacement of H(3 instead of using the high 
ionization line CivA1549. In addition, independent methods based on the well established 
oc <7* and oc M bu i g( , relations for active and non-active galaxies (Gebhardt et al. 2000; 
Ferrarese et al. 2001; Laor 2001; McLure & Dunlop 2001,2002; Tremaine et al. 2002) result in 
consistent black hole mass estimates (see Barth 2003; Vestergaard 2004 for a more detailed 
discussion). 

In the following we present the results of SMBH mass estimates for luminous quasars 
at redshifts 3.3 < z < 5.1 that we observed in the near infrared to study the evolution of the 
Fell/Mgll emission line ratio (Dietrich et al. 2002a,2003a). We analyzed the broad emission 
line profiles of MgllA2798 and Rf3. For most of the quasars we also have observed or have 
access to the ultraviolet spectral range that contains the strong CivA1549 emission line. 
Hence, we could use almost always at least two emission line profiles to determine M bh and 
for the first time to compare the obtained SMBH masses. 

In the following we assume H Q = 65 kms _1 Mpc _1 , Vl M = 0.3, = 0.0. Introducing 
Ha = 0.7 (Netterfield et al. 2002) instead of = 0.0, the luminosities at the highest redshifts 
would be ~ 10 % larger. Hence, a slight change in luminosity will only have a minor impact 
on the estimated black hole mass. 



2. The High-Redshift Quasar Sample 

We have observed a sample of 15 luminous high-redshift quasars with z ~ 3.3 up to 
z ~ 5.1 (Table 1). The spectra were obtained using infrared and optical spectrographs of 
the Cerro Tololo Inter- American Observatory at Chile (OSIRIS, 4 m), the ESO observatories 
at La Silla (Soft, 3.5 m) and Paranal (ISAAC, 8.2 m) at Chile, the Calar Alto Observatory 
(Twin, 3.5 m) at Spain, and the W.M. Keck Observatory (NIRSPEC, 10 m) at Hawai'i (for 
more details on the observations see Dietrich et al. 1999, 2002a, 2003a; Dietrich & Wilhelm- 
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Erkens 2000). In order to study the broad ultraviolet Fell emission and the MgllA2798 
emission line ratio J, H, and K-band spectra were recorded. For the quasars with z ~ 3.5 
we observed in the K-band the rest frame optical wavelength range which covers H/3 and 
[O ill] A4959, 5007. Furthermore, for 10 of the studied quasars we have ultraviolet spectra of 
the CivA1549 emission line (Dietrich et al. 2002b). 



3. Data Analysis 

3.1. Emission Line Profile Measurements 

For each of the high-redshift quasars we measured the emission line profile width of H/3, 
MgllA2798, and CivA1549 if covered. To determine a reliable profile width it is necessary to 
identify and to remove potential contamination by additional emission features. Although 
the MgllA2798 emission line doublet appears to be relatively isolated in quasar spectra, it is 
embedded in the superposition of many thousands of discrete ultraviolet Fe II emission lines 
that form broad emission blends (e.g., Wills, Netzer, & Wills 1985). The careful correction 
for ultraviolet Fe emission yields a more solid measurement of the Mgll profile width than 
assuming a linear pseudo-continuum beneath the Mgll emission line. The contamination of 
B./3 by optical Fell emission is less severe (e.g., Boroson & Green 1992). CivA1549 suffers 
only minor contamination of weak emission lines, like NlV"]A1486 in the blue wing and from 
Fell emission beneath the Civ line profile (e.g., Vestergaard & Wilkes 2001). To account 
for these emission line profile contaminations we analyzed the high redshift quasar spectra 
using a multi-component fit approach (Dietrich et al. 2002a,2003a). Following Wills, Netzer, 
& Wills (1985) we used the following four components to reconstruct each quasar spectrum, 
(i) a power-law continuum (F„ ~ u a ), (ii) Balmer continuum emission (Grandi 1982; Storey 
& Hummer 1995), (iii) Fell emission templates for the ultraviolet and optical wavelength 
range (Vestergaard & Wilkes 2001; Boroson & Green 1992), and (iv) and Gaussian profiles 
to account for individual broad emission lines, such as H/5, Mgll, and Civ. The power-law 
continuum fit to the individual quasar spectra was iteratively determined taking into account 
the measured strength of the Fe II emission and Balmer continuum emission derived by the 
fitting process. The emission line profile widths, expressed as full width at half maximum 
(FWHM) which we measured for H/3, Mgll, and Civ if are accessible, are given in Table 2, 
together with the spectral slope a of the power-law continuum. Since the spectral resolution 
of the infrared spectra is of the order of ~ 200 to 500 kms -1 (Table 2), a correction of the 
measured profile width with respect to the instrumental resolution is not necessary. 
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3.2. Estimating Black Hole Masses 

The empirically justified radius-luminosity relation, Rblr oc L 13 (e.g., (3 = 0.7, Kaspi et 
al. 2000) and the assumption of virial dominated motion of the broad line region (BLR) gas 
(Peterson & Wandel 1999; Onken & Peterson 2002) has led to equations which connect the 
black hole mass M bh with observationally accessible parameters like luminosity and emission 
line profile width. We applied the equations given by Kaspi et al. (2000), McLure & Jarvis 
(2002), Vestergaard (2002), and Warner et al. (2003) to estimate the black hole mass for 
the high-redshift quasars we observed. The explicit form of the equations we employed, are 
presented in the appendix for R/3, MgllA2798, and CivA1549. 

It should be kept in mind that although CivA1549 is a high-ionization emission line, 
in contrast to R(3 and MgIlA2798 as typical low-ionization lines, there are good reasons 
to assume that the broad line Civ emitting gas motion also is gravitationally dominated 
(Peterson & Wandel 2000). This assumption is based on variability response, the emission 
line profile shape, and the relation of the profile widths of H(3 and C IVA1459 combined with 
variability time lags (Onken & Peterson 2002). Furthermore, it has been shown that using 
emission line profiles of H(3 and of C IVA1549 to estimate the SMBH mass of quasars yields 
consistent results on average (e.g., Vestergaard 2002; Warner et al. 2003). 

4. Results 

4.1. Super Massive Black Hole Mass Estimates 

We applied the equations (Al) to (A5) given in the appendix together with the emission 
line profile widths of H/3, Mgll, and Civ and the spectral index a of the continuum slope 
(Table 2), to calculate the mass of the SMBH for each of the high-redshift quasars. To 
estimate the uncertainty of the determined SMBH masses the errors of the emission line 
profile width measurements and of the continuum luminosity have been propagated. In 
Table 3 the resulting black hole masses for each emission line estimator are presented. In 
Figure 1 the black hole mass is displayed for each emission line as a function of redshift. 
The luminous quasars under study have black hole masses in the range of M&/, ~2x 10 8 M Q 
up to M bh ~ 4 x 1O 1O M . The SMBH mass based on the H/5 profile width are consistent 
within a factor of ~ 1.9±0.1, with lower mass estimates provide by the relation suggested by 
McLure & Jarvis (2002) compared to those of Kaspi et al. (2000). This difference is mainly 
caused by the different slopes of the applied R L relation. Assuming a slope of (3 = 0.6 
instead of (3 = 0.7 will result in ~ 2 times lower estimates. Utilizing the Civ emission 
line profile width, the resulting SMBH mass estimates, using equations A3 and A4 as given 
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in the appendix, are consistent within ~ 15 % , as expected since the same slope is assumed 
for the R L relation (eqs. A3 and A4). The scatter of the two SMBH estimates based on 
C IV are due to the use of the individual spectral index a for each quasar that we determined 
instead of using a uniform continuum slope of a — —0.4 for example. Furthermore, to use 
the continuum luminosity L(1450 A) for each quasar an additional a-dependent factor is 
introduced. In spite of these modifications the resulting black hole mass estimates, using 
the C IV profile width, agree well with each other. Close inspection of Figure 3 obviously 
shows that the M bh estimates based on Mgll are systematically lower than the ones based 
on H/3 and CivA1549. On average the MgllA2798 emission line profile indicates a ~ 5x 
lower SMBH mass than the H/3 and CivA1549 line profile. For each of the high-redshift 
quasars we calculated a mean SMBH mass based on H/3, Mg II, and C IV if available (Table 
3). Although the SMBH mass estimates based on Mgll are lower than those given by H/3 
and Civ, we do not apply possible corrections (see Sect. 4.2). 



4.2. Size of the Mgll Emission Region 

A basic assumption of the R-L relation involving the MgllA2798 emission line is that 
the emission regions of Mgll and H/3 are nearly at the same distance from the central 
continuum source. While the comparison of the profile width of both emission lines indicate 
that Rblr{H (3) is nearly identical with i?BLj?(MgII) (McLure & Jarvis 2002), Corbett et 
al. (2003) find that Mg II has generally a broader profile than H/3. In Figure 2 we compare 
the FWHM of the line profiles of Mgll and H/3 for the seven quasars at z ~ 3.5. With the 
exception of Q 2050-359 and Q 0105-2634 (similar width) and Q 2348-4025 (larger Mgll 
profile than H/3) the H/3 profile is generally slightly broader than the Mgll profile. The 
average ratio of the FWHM of the H/3 to MgllA2798 line profile amounts to 1.2 ± 0.5 based 
on the seven z ~ 3.5 quasars we studied. 

The most direct measure of the BLR radius for Mg II is given by reverberation mapping 
experiments. So far, almost all monitoring campaigns dedicated to AGNs have been focused 
in the optical on the H/3 region and in the ultraviolet on the Lya - CivA1549 - Cm]A1909 
wavelength range. However, for NGC5548 (Clavel et al. 1991; Peterson et al. 1991; Dietrich 
& Kollatschny 1995) and NGC3783 (Reichert et al. 1994; Stirpe et al. 1994) MgllA2798 
and H/3 were measured. NGC3783 showed only minor Mgll emission line flux variations 
during the seven month campaign in 1992. The derived emission line flux variation delays 
for Mgll and H/3 are of the same order of r ~ 6 — 8 days (Reichert et al. 1994; Stirpe et 
al. 1994). In contrast, NGC5548 underwent quite strong variations for H/3 and Mgll during 
the first year of its monitoring in 1989. Dietrich & Kollatschny (1995) used the NGC5548 



-7- 



observations taken with IUE in the ultraviolet (Clavel et al. 1991) and from the ground-based 
monitoring campaign (Peterson et al. 1991) to estimate a typical radius of the MgllA2798 and 
H/5 emission line regions. They determined time delays based on the centroid (r c ), as well 
as on the peak (r p ), of the interpolated cross-correlations function (ICCF). The time delays 
based on the ICCF centroid amount to r c (H/3) = 19.7 ± 2.4 d and r c (MgII) = 46.5 ± 2.3 d, 
while the peak of the ICCF yield delays of r p (H/3) = 19.5 ± 2.4 d and r p (MgII) = 39.5 ± 3.3 d 
(Dietrich & Kollatschny 1995). Assuming that these delays represent typical distances of the 
line emitting gas from the central variable continuum source, in the case of NGC 5548 the 
Mg II gas is located ~ 2 times more distant than the H/3 emitting region. In a very simplified 
model of gravitationally bound Keplerian orbits it can be expected that the typical velocity 
of Mgll is about \/2 times narrower than for However, this estimate of R B L R (Mgll) 
is based on only one Seyfert galaxy. Although this result is quite consistent with the ratio 
of the profile widths of H/3 to Mgll that we find for the high-redshift quasars (Figure 2 
and Table 2), it is necessary to measure R B iR(M.gu) for more Seyfert galaxies and quasars. 
Since a larger distance of the Mgll BLR emitting region than for H/3 will introduce a factor 
of b = R BLR (Mgll) j R BLR (R(3), the estimated SMBH mass will be slightly underestimated 
assuming R BLR (Mgll) ~ R BLR (K{3). 



4.3. Bolometric Luminosity and Eddington Ratio 

To compare the SMBH mass with the bolometric luminosity, L bo i, and to determine 
the Eddington luminosity ratio, L bo i/L edd , we calculated L bo i. To transform measured rest 
frame luminosities \L\ into L boi scaling factors are provided by e.g. Elvis et al. (1994) and 
Kaspi et al. (2000). Most recently, Vestergaard (2004) and Warner et al. (2004) present 
revised scaling factors for optical and ultraviolet continuum luminosities employing results 
of multi-frequency studies on the AGN spectral energy distributions from X-ray energies to 
the mid-infrared wavelength range. They found consistent scaling factors which are slightly 
smaller (~ 20%) than previously assumed. We used the relation L bo \ = 4.36 x AL^(1450A) 
to calculated the bolometric luminosity (Warner et al. 2004) for each of the high-redshift 
quasars under study (Table 1). In Figure 3 we plot the average SMBH masses as a function 
of Lboi- At a given L bo i the determined M b h of the high-redshift quasars cover a range of up to 
one order of magnitude in mass, similar to comparable studies for quasars at lower redshift 
(e.g., Woo & Urry 2003; Vestergaard 2004; Warner et al. 2004). 

For a fully ionized pure hydrogen gas, i.e. \i e = 1, the Eddington luminosity is given as 
(Rees 1984; Peterson 1997) 
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L edd = 1.26 x 10 38 (^) ^ erg s" 1 (1) 

Assuming a more realistic gas chemical composition containing hydrogen, helium, and 
some heavier elements will result in an additional factor /j, e = 1.15 accounting for the atomic 
weight per electron (Haiman & Loeb 2001). This will yield a 15% higher Eddington lumi- 
nosity. In the following we will assume pure hydrogen gas because /i ^ 1 has only minor 
impact on L e dd- 

The Eddington ratio L bo i/L e dd is determined for each quasar individually. The mean 
black hole masses are used to calculate the Eddington luminosity, L e dd, for spherical accre- 
tion. Together with the bolometric luminosity, the Eddington ratio, Lb i/L e dd, is calculated 
(Table 3). The determined L bo i/ L e dd covers a range of ~ 0.5 to ~ 7 which is in good agree- 
ment with results of studies of quasars sample that are dominated by low and intermediate 
redshift quasars (e.g., Woo & Urry 2003; Vestergaard 2004; Warner et al. 2004). Hence, 
there is no difference regarding the Eddington ratio L bo i/L e dd for quasars at high-redshifts 
compared to those at lower redshift. 



5. Discussion 

5.1. Comparison of M bh based on H(3, Mgll, and C IV 

In Table 3 the SMBH mass estimates based on the three emission lines H/5, MgllA2798, 
and CivA1549 are listed for comparison. Similar to previous studies there appears to be 
an approximate upper limit for at ~ 10 10 M Q (e.g., Corbett et al. 2003; Netzer 2003; 
Shields et al. 2003; Vestergaard 2004; Woo & Urry 2003; Warner et al.2003). Only Q0105- 
2634 yields an estimate of M bh = (4.1 ± 1.2) x 1O 1O M (eq. Al) and (2.1 ± 0.4) x 1O 1O M 
(eq. A2), respectively, based on the H/5 emission line width. However, the H/5 emission profile 
of this particular quasar has a rather low signal-to-noise ratio. Furthermore, the estimate 
M bh = (4.1 ± 1.2) x 1O 1O M based on the R-L relation with a slope of (3 = 0.7. Assuming 
a flatter slope f3 of the R-L relation, as suggested by Peterson et al. (2000) and McLure & 
Jarvis (2002) the mass estimate is reduced by a factor of ~ 2 and in this case within the 
uncertainties in better agreement with an empirical upper limit of ~ 10 10 M Q which 

also has been suggested by Umemura (2003) and Rocca-Volmerange et al. (2004) who studied 
early phases of galaxy formation and related black hole growth. 

For two quasars at redshift z ~ 3.5 we detected quite prominent [O ill] A4959, 5007 
emission. It has been suggested that the emission line profile width of [O ill] A4959, 5007 
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can be used as a surrogate of the stellar velocity dispersion <r* (Nelson 2000). Applying the 
relation that connects black hole mass with stellar velocity dispersion (Tremaine et al. 2002) 
the use of [O ill] A4959, 5007 as a SMBH mass estimator has been investigated by Shields et 
al. (2003). Using the emission line profile width of [O in] A5007 the estimated the mass of the 
black holes amount to 1.4 ±0.3 x 1O 9 M (Q 0256-0000) and 1.5 ±0.3 x 1O 9 M (Q 0302-0019) 
which are consistent with the masses provided by Civ, Mgll, and H(3 (Table 3). 

While the black hole mass estimates based on CivA1549 and H/3 are quite consistent, 
the estimate provided by MgllA2798 indicates lower masses (Table 3). Although R - 
L relations based on statistical results found for quasar samples introduce uncertainties of 
a factor of at least ~ 3 in the case of individual quasars (e.g., Vestergaard 2002) there are 
some additional possibilities which may cause lower M bh estimates based on MgllA2798. 

First, the slope (3 of the R oc L 13 relationship is still not well defined, varying between 
~ 0.5 to 0.7 (e.g., Kaspi et al. 2000; Peterson et al. 2000; McLure & Jarvis 2002; Vestergaard 
2002; Netzer 2003). A flatter slope in the R — L relation will result in less massive black holes. 
While in the original analysis by Kaspi et al. (2000) it is found that (3 = 0.70 ± 0.03, there 
are indications that the slope is closer to (3 ~ 0.62 ± 0.02 (Peterson et al. 2000). McLure & 
Jarvis (2002) derived (3 = (0.47 ± 0.05) for the R - L relation using Mgll. They determined 
this flatter slope under the assumption that Mg II and B.f3 originate at nearly the same BLR 
radius. Hence, they used B./3 reverberation results for Seyfert galaxies and quasars to analyze 
the correlation of rBi,i?(Mgll) with the luminosity ALa(3000). However, assuming a slightly 
flatter slope of (3 = 0.6 instead of (3 = 0.7 will decrease the estimated black hole mass based 
on H/5 by a factor of ~ 2. Using furthermore a slightly steeper slope for Mgll (e.g., (3 = 0.6 
instead of (3 = 0.47) will provide on average ~ 2.4 ± 0.2 times higher SMBH estimates. 
A second aspect as discussed in Sect. 4.2 is the size of the Mgll emitting region. If Mgll 
originates at slightly larger distances to the center of the continuum emission than H/3, the 
estimated mass can be increased by an additional factor. Currently, we suppose that the 
different SMBH estimates provided by Civ and H/3 compared to those based on Mgll arc 
predominantly caused by the different slopes of the R-L relations. Furthermore, it is still 
an open question whether the R-L relation can be extrapolated to quasar luminosities since 
the relation is derived for active galactic nuclei with L hoi < 10 46 ergs -1 (e.g., Corbett et 
al. 2003). In the case of quasars this relation is extended to ~ 10 2 to ~ 10 3 times more 
luminous objects. Yet, in spite of the uncertainties regarding the use of the H/3, Mgll, and 
C IV emission line profiles to determine the black hole mass, we emphasize that all three 
emission lines indicate consistently the presence of super-massive black holes with masses of 
10 9 < M bh < 1O 1O M at redshifts up to z=5.1. 
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5.2. Implications of Black Hole Growth Time Scales 



The presence of SMBHs at high-redshifts and the close relation of quasar activity and 
galaxy formation provide valuable constraints on the epoch when these SMBHs had to start 
to grow in their forming host galaxies. The time span that is necessary to build-up a 
~ lO 9 tolO lo M SMBH can be described by an e-folding time for the growth of a single 
seed black hole (e.g., Haiman & Loeb 2001). We use the following equation to estimate the 
e-folding time scale and hence the epoch when the seed black hole started to grow (Salucci 
et al. 1999) 



with r = t obs — t as time elapsed from the initial time, t , to the observed time, t t, s , M^ ed 
as seed black hole mass, 77 = L bo i/ L e dd as Eddington ratio, e as efficiency to convert mass to 
energy, and as Eddington time scale, with = <Jtc/ 'AnGrrip = 3.92 x 10 8 yr (Rees 1984). 
The Eddington time t# describes the time necessary to radiate at the Eddington luminosity 
the entire rest mass of an object. In hierarchical models for structure formation the first 
baryonic objects which collapsed had masses of the order of ~ 10 4 — 1O 6 M (e.g., Silk & Rees 
1998; Larson 2000; Shibata & Shapiro 2002; Bromm & Loeb 2003). This provides upper 
limit estimates for M^ ed of black holes. Recent models of early star formation indicate that 
nearly metal free Pop III stars, formed at high-redshifts (z > 20), were predominantly very 
massive with M > 100M Q (e.g., Fryer, Woosley, & Heger 2001; Abel, Bryan, & Norman 
2002; Bromm, Coppi, & Larson 2002). The stellar remnants are expected to be of the order 
of ~ 1OM (e.g., Fryer 1999). However, some early star formation models even indicate the 
possibility of black hole remnants with several times ~ 1O 3 M (Bond, Arnett, & Carr 1984). 
In the following the efficiency e to convert mass to energy is assumed to be e = 0.1. It should 
be noted that the efficiency e can reach values up to e = 0.42 for a maximal rotating Kerr 
black hole. 

We employed equation (2) to derive the time r which is necessary to accumulate the 
determined mean black hole mass (Table 3) for seed black holes with M^ ed = 1OM , 1O 3 M , 
and 1O 5 M , respectively. We find that the time r to build-up a SMBH is between several 
100 Myr and ~ 1 Gyr, i.e., almost always r is smaller than the age of the universe at the 
corresponding redshift. To illustrate this in a proper way we used r to calculate t , i.e. the 
epoch when the SMBH, observed at t obS) started to form. In Figure 4 the time t is displayed 
as a function of the observed redshift z obs . The determined time scales r and hence t indicate 
that most of the SMBH residing in the observed high-redshift quasars started to grow at 
redshifts of z(t ) ~ 5 to 10. As can be seen for low massive black hole seeds, i.e., stellar 




(2) 
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remnants with M^ ed ~ 1OM some black holes had to have started to grow at very early 
epochs and for Q 0105—2634 and SDSS 0388+0021 the required time ris even larger than the 
actual cosmic age at z obs . To avoid too long growing times r, higher M^ ed masses appear to 
be favored. Assuming high mass seed black holes with M^ ed = 10 3 to 10 5 M Q , the SMBH in 
the z ~ 3.5 quasars began to grow at redshifts z(t D ) ~ 4 and earlier. The luminous quasars 
at redshifts z > 4.5 show a smaller scatter for t Q compared to those at z ~ 3.5. Based on our 
results the SMBH residing in the z > 4.5 quasars started at z(t Q ) ~ 5 to 10. However, stellar 
mass black hole seed are still a valid option because a flatter slope of the R L relation for 
H/? with (3 ~ 0.6 (Peterson et al. 2000; McLure & Jarvis 2002) will yield about two times less 
massive black holes. Hence, the Eddington luminosity will be reduced and the Eddington 
ratio will be increased. As a result the growing time scale r will decrease and SMBHs can 
be built-up with stellar mass remnants as black hole seeds as well. 

It is remarkable that most of the super-massive black holes of the studied high-redshift 
quasars are in place at the end of the era of reionization of the universe which is supposed 
to end at z ~ 6 (e.g., Becker at al.2001; Fan et al. 2002). Particularly, for the quasars at 
redshifts z > 4.5 we want to emphasize that the epoch when these black holes started to grow 
coincides with the beginning of first intense star formation that marks early phases of the 
host galaxy formation. This epoch of major star formation also is indicated by the chemical 
enrichment history of the quasar gas. To achieve the estimated metallicities of several times 
solar (e.g., Dietrich et al. 2003b, c) a major and intense star formation episode had to start at 
z ~ 6 to 8. Further circumstantial evidence for a vigorous star formation phase at those high- 
redshifts is given by the lack of evolution of the iron/a element ratio compared to results 
for local quasars (e.g., Barth et al. 2003; Dietrich et al. 2002,2003a; Freudling et al. 2003; 
Iwamuro et al. 2002; Maiolino et al. 2003; Thompson et al. 1999). 

The presence of super-massive black holes in luminous high-redshift quasars, the esti- 
mated time scale required to build-up black holes of several M bh ~ 1O 9 M , the gas metallicity 
of at least solar up to several times solar, and the lack of evolution of the relative iron/a 
element abundance ratio provide evidence that those quasars reside in already mature host 
galaxies, at least with an evolved massive central component (Hamann et al.2004). While 
the massive galaxy is forming, i.e., major and vigorous star formation activity occurs, the 
super-massive black hole grows. It takes about ~ 0.5 Gyr to accumulate a SMBH with 
several times M bh ~ 1O 8 M to M bh ~ 1O 9 M that is necessary to power quasar activity. 
It takes approximately this same amount of time for a stellar population to enrich the gas 
around quasars to the observed solar or super-solar levels (e.g., Dietrich et al. 2003b, c; Mat- 
teucci & Padovani 1993; Matteucci & Recchi 2001). Therefore, by the time most quasars 
become visible, even at the highest redshifts, a substantially evolved stellar population is 
harboring the active nucleus. This conclusion is supported by some theoretical models of 
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joint SMBH - host galaxy evolution (Archibald et al. 2002, Kawakatu et al. 2003). Direct 
imaging studies of low redshift quasars show clearly that the host galaxies are characterized 
by at least moderately old, stellar populations on >kpc scale (Nolan et al. 2001; Dunlop et 
al. 2003). Although there is less imaging data at higher redshifts at least some quasars still 
have substantial hosts (Kukula et al. 2001). 

It can be assumed that early phases of galaxy formation that are accompanied by intense 
and vigorous star formation also mark the era when super-massive black holes are growing 
from massive seed remnants (Hamann et al. 2004). A joint era of SMBH growth and early 
phases of galaxy formation can be understood in the context of recently suggested galaxy 
formation models (e.g., Archibald et al. 2002; Di Matteo et al. 2003; Granato et al. 2001,2004). 
These models further suggest that the M bh oc relation evolves over cosmic time scales and 
is finally established in the local universe (Di Matteo et al. 2003). A less steep M bh oc a* 
relation also has been suggested by Shields et al. (2003). This may provide an explanation to 
avoid huge host galaxy masses (~ 10 13 M Q ) and large stellar velocity dispersions (cr* ~ 700 
kms -1 ) which have never been observed (e.g., Netzer 2003). 



6. Conclusion 

We have analyzed a sample of 15 luminous high-redshift quasars (3.3 < z < 5.1) to 
estimate the mass of super-massive black holes. For the first time results for M bh based on 
three different emission lines, Civ, Mgll, and H/3, are compared for high-redshift quasars. 

We determine black hole masses in the range of M bh ~ 2x 1O 8 M up to M bh ~ 4x 1O 1O M 
Similar to prior studies there appears to be an approximate upper limit for M bh at ~ 1O 1O M . 
The black hole mass estimates based on C IV and H/3 are quite consistent while Mg II indicates 
by a factor of ~ 5 lower SMBH masses. In spite of the uncertainties regarding the use of 
the H/3, Mgll, and Civ emission line profiles to determine the black hole mass, all three 
emission lines indicate consistently the presence of super-massive black holes with masses on 
average ~ 6 x 1O 9 M at redshifts up to z=5.1. 

We determine an Eddington luminosity ratio of L bo i/L e dd in the range of ~ 0.5 to ~ 7. 
This is consistent with results of studies of quasars sample that are dominated by low and 
intermediate redshift quasars. This indicates that high-redshift quasars are not different 
regarding the Eddington ratio L bo i/L edd compared to those at lower redshift. 

To date the epoch when the SMBHs which drive high-redshift quasar activity should 
have started to grow, we assume logarithmic growth by spherical accretion with a mass to 
energy conversion efficiency of e = 0.1. Generally, the growth-time r is of the order of several 
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~ 100 Myr Assuming higher mass seed black holes {M^ ed = 10 3 tol0 5 M Q ) the SMBH in 
the z ~ 3.5 quasars begun to grow at redshifts z(t ) ~ 4 and earlier while the quasars with 
z > 4.5 started at z(t Q ) ~ 5 to 10. 

The presence of super-massive black holes in luminous high-redshift quasars, the esti- 
mated time scale required to build-up black holes of several times M bh ~ 1O 9 M , the several 
times solar gas metallicity and the lack of evolution of the relative iron/a element abundance 
ratio provide evidence that luminous quasars at high-redshift reside in already mature host 
galaxies, at least with an substantially evolved massive central stellar component. This re- 
sult indicates that early phases of massive galaxy formation that are accompanied by intense 
and vigorous star formation also mark the era when super-massive black holes are growing 
from massive seed remnants. The estimated black hole growth time scale suggest a delay 
of about ~ 0.5 Gyr of the quasar activity with respect to preceding intense star formation 
which is consistent with the observed high gas metallicity. 

We dedicate this Paper to the memory of Richard J. Elston, who left us far too early 
on January 26, 2004. 

MD and FH acknowledge support from NASA grant NAG 5-3234 and NSF grant AST- 
99-84040 (University of Florida). 

A. Appendix 

The equations we applied, given by Kaspi et al. (2000), McLure & Jarvis (2002), Vester- 
gaard (2002), and Warner et al. (2003) have been modified for the use of the continuum 
luminosity L\ that we measured at A = 1450 A for each of the high-redshift quasars of this 
study. We explicitly present the dependence on the spectral index a which transforms the 
continuum luminosity measured at different wavelengths. The uniform use of L A (1450A) 
does not introduce an additional uncertainty, because instead of assuming a fixed continuum 
slope a to rescale L\ as in prior studies we used the measured continuum slope a for each 
quasar which we determined for the power-law continuum fit (F u oc u a , Table 2) in the 
multi-component analysis of the quasar spectra (Dietrich et al. 2002a,2003a). The contin- 
uum slope is measured for a wavelength range of AA ~ 1100 to 5500 A (z = 3.5 quasars) 
and AA ~ 1100 to 4300 A (z = 4.5 quasars). 

For the H/3 emission line profile we determined Mbh applying the equation given by 
Kaspi et al. (2000) (eq. Al) and presented by McLure & Jarvis (2002) (eq. A2). The main 
difference is given by the slope of the radius - luminosity (R - L) relation. 
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M bh = 4.82 x 10» x (mf (l+ "> x (W)" x (FWHM( W 
V5100/ V 10 44 erg/s / V 10 3 km/s / 



The black hole mass is estimated based on the CivA1549 line profile width applying 
the following equations (Vestergaard 2002, eq. A3; Warner et al. 2003, eq. A4). Warner et 
al. (2003) assumed that the radius of the C IV BLR is half the size of the H/5 emitting region 
which is well justified by variability studies (Peterson 1993). 



^q* * (-> 

fi /1450\°- 7 ( 1+a ) /AL A (1450A)\ - 7 /FWHM(CIV)\ 2 „, , 4 „, 

M bh = 2.41 x 10 6 x x — ^ t-L x — -— — '- M G A4 

bh V5100/ V 10 44 erg/s / V 10 3 km/s / o \ } 

To estimate based on the MgllA2798 emission line profile width we employed the 
relation given by McLure & Jarvis (2002). 



V3000/ V 10 44 erg/s / V 10 3 km/s / o \ j 
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Fig. 1. — Estimates SMBH masses for the high-redshift quasars as a function of redshift 
using eq. (Al) to (A5). The triangles indicate based on H/5 (A - Kaspi et al. 2000; 
▲ - McLure & Jarvis 2002), the squared symbols show based on Civ (□ - Warner 
et al. 2003; ■ - Vestergaard 2002), and the filled diamonds display M bh based on Mgll (♦ 
McLure & Jarvis 2002). 

Fig. 2.— Comparison of the line profile width, FWHM, of MgllA2798 and H/3 for the 
high-redshift quasars at z ~ 3.5. The dotted line indicates a perfect one-to-one relation. 

Fig. 3. — Mean SMBH masses for the high-redshift quasars as a function of bolometric 
luminosity, Lw, are shown. The displayed errors are taking into account the measurement 
uncertainties and the distribution of the individual SMBH mass estimates in quadrature 
(Table 3). 

Fig. 4. — The time t Q = r obs — r is plotted as a function of redshift for several seed black 
hole masses to illustrate when the super-massive black holes had had to start to form. The 
three panels display results for different seed black hole masses with M^ ed = 10M Q (left), 
M^ ed = 1O 3 M (center), and M^ ed = 1O 5 M (right). The redshift which corresponds to t 
is labeled at the right as z(t a ) (H G = 65 kms -1 Mpc" 1 , VL M = 0.3, = 0.0). 
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Table 1. The high-redshift quasars have been observed at the observatories and dates as 
given in (2) and (3). In addition, the achieved spectral resolution for the emission line 

profiles of H/3, Mgll, and Civ is listed in unit of kms -1 (columns 4 to 6). The redshift of 
the quasars is given in (7). The rest frame continuum luminosity AL A (1450A) and the 
bolometric luminosity, L^, are listed in (8) and (9). 



quasar 


site 


date 




resol. 




redshift 


log AL A (1450A) 


iogL bol 








CIV 


Mgll 


H/3 
















[kms- 1 ] 






[ergs -1 ] 


[ergs^ 1 ] 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


BRI 0019-1522 


CTIO 


09-00 


300 


230 




4.53 


46.92 


47.56 


BR 0103+0032 


CTIO 


09-00 


300 


230 




4.44 


46.99 


47.63 


Q 0103-260 


VLT 


12-98 


750 






3.38 


46.83 


47.47 




NTT 


10-99 




500 


430 








Q 0105-2634 


NTT 


10-99 




500 


430 


3.49 


47.46 


48.10 


PSS J0248+1802 


CTIO 


09-00 


300 


230 




4.44 


47.24 


47.88 


Q 0256-0000 


NTT 


10-99 


300 


500 


430 


3.38 


47.14 


47.78 


Q 0302-0019 


NTT 


09-99 


300 


500 


430 


3.29 


47.24 


47.88 


SDSS 0338+0021 


VLT 


12-98 


750 






5.00 


46.61 


47.25 


PC 1158+4635 


KECK 


05-00 


300 


520 




4.73 


47.39 


48.03 


SDSS 1204-0021 


VLT 


04-01 




525 




5.09 


46.79 


47.42 


Q 2050-259 


CTIO 


09-00 




230 


230 


3.51 


47.44 


48.08 


PKS 2126-158 


CA 


08-93 


300 






3.27 


47.61 


48.25 




CTIO 


09-00 




230 


230 








Q 2227-3928 


NTT 


10-99 




500 


430 


3.44 


47.07 


47.71 


BRI 2237-0607 


CTIO 


09-00 


300 


230 




4.56 


47.27 


47.91 


Q 2348-4025 


NTT 


10-99 




500 


430 


3.31 


47.27 


47.91 
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Table 2. Measurements of the emission line profile widths (FWHM) for H/3, MgllA2798, 
and CivA1549 are given for the studied high-redshift quasar sample. The slope a of the 
power-law continuum fit (F v oc u a ) given in column (5). 



quasar 


FWHM(H/3) 


FWHM(Mgn) 
[kms- 1 ] 


FWHM(Civ) 


a 


(1) 


(2) 


(3) 


(4) 


(5) 


BRI 0019-1522 




1993 ± 353 


4078 ± 347 


+0.10 ±0.09 


BR 0103+0032 




1869 ± 268 


2626 ±118 


-0.46 ± 0.03 


Q 0103-260 


4384 ± 432 


2725 ± 193 


3257 ± 350 


-1.33 ±0.06 


Q 0105-2634 


6537 ± 650 


6377 ± 570 




-0.69 ±0.06 


PSS J0248+1802 




3811 ± 279 


5620 ± 224 


-0.48 ±0.08 


Q 0256-0000 


4247 ± 277 


3055 ± 182 


3784 ± 163 


-0.41 ±0.05 


Q 0302-0019 


3257 + 191 


2457 ± 171 


2102 ± 108 


-0.21 ±0.03 


SDSS 0338+0021 






4500 ± 1200 


-0.01 ±0.20 


PC 1158+4635 




4559 ± 257 


4833 ± 300 


-0.37 ±0.05 


SDSS 1204-0021 




6090 ± 380 




-0.43 ±0.06 


Q 2050-259 


4589 ± 247 


4335 ± 375 




-0.15 ±0.18 


PKS 2126-158 


2707 ± 320 




4991 ± 176 


-1.09 ±0.09 


Q 2227-3928 


4812 ± 326 


2890+ 118 




-0.85 ±0.05 


BRI 2237-0607 




2517 + 214 


3550 ± 269 


-0.55 ±0.02 


Q 2348-4025 


1061 ± 142 


4315 ± 193 




-0.53 ±0.07 
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Table 3. The SMBH masses estimates for the studied high-redshift quasar sample. The 
derived black hole masses based on H/5 using Kaspi et al. (2000, K00) and McLure & Jarvis 
(2002, MJ02) are given in (2) and (3). The SMBH estimates employing CivA1549 are listed 
in column (4) and (5) ( Vestergaard 2002, MV02; Warner et al2003, CW03) and using 

MgllA2798 in (6) (McLure & Jarvis 2002, MJ02). In column (7) the mean super-massive 
black hole mass is listed and the resulting Eddington ratio Lf, i/L e dd is given in column (8). 



quasar 


H/3 (K00) 


H/3 (MJ02) 


C iv (MV02) 


C iv (CW03) 


Mgn (MJ02) 


M Wi 


Lfcoi/Ledd 








M bh [1O 9 M ] 






[1O 9 M ] 




(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


BRI 0019-1522 






3.1 ± 0.6 


1.7 + 0.3 


0.22 ± 0.08 


1.7 + 0.6 


1.73 


BR 0103+0032 






1.4 + 0.4 


1.3 + 0.4 


0.25 ± 0.07 


1.0 + 0.3 


3.48 


Q 0103-260 


11.8 + 2.0 


6.2 + 1.3 


1.6 + 0.4 


3.3 + 0.8 


0.60 ± 0.10 


4.7+1.5 


0.50 


Q 0105-2634 


41.4 + 12.2 


20.6 + 4.5 






5.2 + 1.0 


22.4 + 8.0 


0.45 


PSS J0248+1802 






9.6+ 1.3 


8.9 + 1.2 


1.4 + 0.2 


6.6+1.9 


0.91 


Q 0256-0000 


8.1 + 2.0 


4.5 + 0.6 


3.7+ 1.4 


6.9 + 0.4 


0.77 + 0.10 


4.8 + 0.9 


0.99 


Q 0302-0019 


4.7 + 1.3 


2.6 + 0.3 


1.4 + 0.2 


1.0 + 0.1 


0.52 ± 0.08 


2.0 + 0.6 


2.95 


SDSS 0338+0021 






2.3+ 1.2 


2.7 + 1.5 




2.5 ± 0.7 


0.56 


PC 1158+4635 






9.0+ 1.5 


7.6 + 1.2 


2.2 ± 0.3 


6.3+ 1.5 


1.35 


SDSS 1204-0021 










2.1 ± 0.3 


2.1 ± 0.3 


1.00 


Q 2050-259 


12.2 ± 3.9 


6.5 + 0.9 






2.0 + 0.4 


6.9 + 2.3 


1.38 


PKS 2126-158 


12.8 + 3.1 


5.9 + 1.5 


13.1 ± 1.7 


21.7 + 2.7 




13.4 + 2.4 


1.01 


Q 2227-3928 


13.8 + 2.7 


7.3 + 1.2 






0.74 ± 0.08 


7.3+ 1.5 


0.56 


BRI 2237-0607 






4.0 + 0.7 


4.0 + 0.7 


0.63 + 0.12 


2.9 + 0.8 


2.25 


Q 2348-4025 


0.70 + 0.30 


0.37 + 0.10 






1.8 + 0.2 


1.0 + 0.3 


6.67 



